Abstract. We revisit the nature of the far infrared 
(FIR) /Radio correlation by means of the most recent mod- 
els for star forming galaxies, focusing in particular on the 
case of obscured starbursts. We model the IR emission 
with our population synthesis code, GRASIL (Silva et 
al. 1998). As for the radio emission, we revisit the sim- 
ple model of Condon & Yin (1990). We find that a tight 
FIR/Radio correlation is natural when the synchrotron 
mechanism dominates over the inverse Compton, and the 
electrons cooling time is shorter than the fading time of 
(-v^ | the supernova (SN) rate. Observations indicate that both 

■ these conditions are met in star forming galaxies, from 
| normal spirals to obscured starbursts. However, since the 

■ radio non thermal (NT) emission is delayed, deviations 
£h \ are expected both in the early phases of a starburst, when 

, ^ ■ the radio thermal component dominates, and in the post- 
starburst phase, when the bulk of the NT component orig- 
inates from less massive stars. We show that this delay 
I | allows the analysis of obscured starbursts with a time res- 

■ olution of a few tens of Myrs, unreachable with other star 
Q\ j formation (SF) indicators. We suggest to complement the 
O^l ■ analysis of the deviations from the FIR/Radio correlation 

with the radio slope (q-Radio slope diagram) to obtain 
, characteristic parameters of the burst, e.g. its intensity, 
| age and fading time scale. The analysis of a sample of com- 
f^*> pact ULIRGs shows that they are intense but transient 
"^h | starbursts, to which one should not apply usual SF indi- 
O . cators devised for constant SF rates. We also discuss the 
possibility of using the q-radio slope diagram to asses the 
presence of obscured AGN. A firm prediction of the mod- 
^ ' els is an apparent radio excess during the post-starburst 
C$ • phase, which seems to be typical of a class of star forming 
galaxies in rich cluster cores. Finally we discuss how devi- 
ations from the correlation, due to the evolutionary status 
/\ ' of the starburst, affect the technique of photometric red- 
^ . shift determination widely used for the high-z sources. 
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1. Introduction 

In recent years the study of starburst galaxies has become 
a very popular subject because of its intimate connection 
with the global star formation history of the Universe. On 
one side high redshift observations in the optical bands 
probe rest frame spectral regions that are highly affected 
by even tiny amounts of ongoing star formation and dust 
extinction. On the other, theoretical models following the 
paradigm of the hierarchical clustering scenario predict 
that merging induced star formation should be highly en- 
hanced in the past. Current estimates of the star formation 
rate (SFR) of the Universe have thus been interpreted on 
the basis of our understanding of local analogous galax- 
ies, in particular through UV continuum and optical line 
emission. However in local starbursts a significant frac- 
tion of the ongoing star formation may be hidden to UV 
and optical estimators. In fact, though starburst galaxies 
were initially selected for the prominence of their optical 
emission lines, it appears that this criterion excludes other 
actively star forming objects and possibly limits our un- 
derstanding to a small phase of their evolution. After the 
IRAS satellite, it became clear that SF is also highly en- 
hanced in very and ultra luminous infrared galaxies that 
are otherwise highly attenuated in the optical. With space 
densities similar to those of quasars (Soifer et al. 1986) 
and total infrared luminosities spanning the range 10 11 - 
10 12 L Q and above 1O 12 L , respectively, Luminous and Ul- 
traluminous Infrared galaxies (LIRGs and ULIRGs) are 
the most luminous objects in the local Universe. Evidence 
of the important role played by dust reprocessing was pro- 
vided by the detection of a diffuse FIR background whose 
high intensity (equal to or higher than that of the opti- 
cal, e.g. Hauser et al. 1998) implies that these galaxies are 
undergoing intense star formation activity (Puget et al. 
1996; Dwek et al. 1998). Furthermore, the advent of the 
Infrared Space Observatory (ISO), in combination with 
the availability of new ground facilities such as SCUBA 
on the JCMT, have discovered numerous high-z galaxies 
with enhanced IR emission (e.g. Elbaz et al. 1999; Smail 
et al. 2000; Barger et al. 2000). 
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Silva et al. (1998) have first introduced the concept 
of age-selective obscuration, to explain the features of the 
observed spectral energy distribution (SED) of star form- 
ing galaxies, from normal spirals to dust obscured star- 
bursts, from the UV to the sub-mm. In this model young 
stars arc supposed to originate within molecular clouds 
and correspondingly their light is attenuated more than 
that of older stars, that already got rid of their parental 
cloud. The UV light in many starbursts is thus dominated 
by the older stars rather than by the younger popula- 
tions. With the same assumptions, Granato et al. (2000) 
reproduced the observable properties of local galaxies (in 
particular the IRAS luminosity function), working within 
the context of structure formation through hierarchical 
clustering, which has successfully confronted a wide range 
of observations on large scale structure and microwave 
background anisotropics. They showed that the concept 
of age-selective obscuration could explain the difference 
between the galactic extinction law and the attenuation 
law observed in starburst galaxies (Calzetti, Kinney, & 
Storchi-Bcrgmann 1994). 

Poggianti, Bressan, & Franceschini (2001) have re- 
cently investigated the optical spectra of very luminous 
infrared galaxies to constrain the recent history of SF and 
the dust extinction characteristic of various stellar popu- 
lations. They have found that the most plausible explana- 
tion for their unusual combination of strong H<5 absorp- 
tion and moderate [Oil] emission is again age-selective 
extinction. Indeed HII regions (wherein the [Oil] emission 
originates) are highly embedded and thus are affected by 
a greater extinction compared to the older stellar popu- 
lations which are responsible for the Balmer absorption. 
Under standard assumptions for the IMF, the SFR derived 
from the fit of the optical spectrum (continuum, absorp- 
tion and emission lines) may account for a small fraction 
of the FIR emission. Moreover, even complemented with 
the information on the FIR flux, the optical-UV spectrum 
is not enough to identify univocal solutions. Further evi- 
dence along this direction is provided by recent observa- 
tions of UV properties of ULIRGs (Goldader et al. 2002). 

These studies underline an intrinsic difficulty of evalu- 
ating the properties of massive starbursts only from their 
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UV, optical and even NIR properties and the natural way 
out from this impasse seems provided by studies at longer 
wavelengths. The capability of FIR and radio spectral re- 
gions to reveal otherwise hidden complex phenomena in 
star forming galaxies, is testified by the existence of a 
"miraculous" correlation between their properties in these 
spectral windows. The FIR/Radio correlation is locally 
well established over a significant range of luminosity, from 
normal spirals to the most extreme ULIRGs and its small 
scatter states the universal proportions with which en- 
ergy is radiated away at IR and radio wavelengths. In 
spite of its obscure nature, its utility appears in several 
aspects, beside being a firm tool within the manifold of 
star formation indicators. For instance the validity of the 
FIR/Radio correlation has been recently confirmed up to 
redshift ~ 1.3 (Garret 2001) and it is widely extrapo- 
lated much beyond, to estimate the redshift of more dis- 
tant objects (e.g. Carilli & Yun 2000). Also, deviations 
from the correlation observed toward the central regions 
of rich cluster of galaxies, where a significant fraction of 
star forming galaxies show a radio excess, are used to trace 
the effect of the hot intracluster medium on their galac- 
tic magnetic field (Gavazzi & Jaffe 1986; Miller & Owen 
2001). 

So far there have been many attempts to explain the 
FIR/Radio correlation but all have soon or late invoked 
a fine tuning of the relevant physical properties, such as 
the intensity of the radiative and magnetic energy den- 
sity (e.g. Lisenfeld, Volk & Xu 1996). In this paper we 
revisit this correlation by combining our spectrophotomet- 
ric code GRASIL, particularly suited to the study of the 
IR properties of dusty galaxies, with a new model of radio 
emission. The latter essentially follows the recipes by Con- 
don (1992), but after a careful assessment of the validity 
of one of its basic assumptions, namely the proportion- 
ality between the non thermal (NT) radio emission and 
the core-collapse supernova (CCSN) rate. This fact ren- 
ders the FIR/Radio correlation so robust and we provide, 
for the first time, a simple explanation of its universality. 
For the same reason we show that deviations are to be 
expected during and soon after the starburst episode, and 
we suggest that they can help in constraining the star for- 
mation history of these galaxies, something that cannot 
be done with optical, NIR and even FIR observations. 

In Sect. H we briefly describe our population synthe- 
sis code for dusty galaxies. Sect. H. is devoted to the new 
model of radio emission. In Sect.^J we describe our cali- 
bration of the NT radio emission model and obtain new 
relations for the SFR against radio emission for the case 
of quiescent galaxies. In Sect. || we analyse infrared and 
radio properties of starburst galaxies. We show that the 
different fading times of FIR and radio emissions may be 
used to reach a time resolution of a few tens of Myr, which 
is impossible resting only on the UV-FIR. In Sect. |^ we 
introduce a new diagnostic tool, the FIR/Radio (q) vs ra- 
dio spectral slope diagram, which potentially allows the 



determination of the evolutionary status of a starburst in 
absence of a good radio spectral coverage. We examine the 
location of an observed sample of compact ULIRGs and 
discuss whether this diagram may also provide a quan- 
titative estimate of the threshold between AGN and star 
formation powered ULIRGs. In Sect, ^we analyse the evo- 
lution in the post-starburst phase and suggest that radio 
excess is actually an indication of the occurrence of this 
phase rather than an environmental effect. In Sect. ^ we 
discuss the impact of these new findings on the determi- 
nation of the redshift of SCUBA sources, a method that 
relies on the FIR/Radio correlation. Sect. |§ is devoted to 
our conclusions. 

2. Infrared Emission 

Infrared emission is calculated with GRASIL, a code de- 
signed to perform population synthesis in presence of dust 
(Silva et el. 1998). In brief, the star formation history, the 
metal enrichment and the current gas fraction are pro- 
vided by a chemical evolution code. Stars and dust are 
distributed either in a disk or in a bulge or both, includ- 
ing high density clumps, the molecular clouds. Young stars 
are supposed to originate within molecular clouds and to 
leave them in a characteristic time scale t esc . The volume 
emissivity is computed by considering the light of young 
stars, absorbed by the molecular clouds, and the light of 
older populations, both propagated through the diffuse 
dust component. The code has been thoroughly tested 
against observations (e.g. Silva et al. 1998; Granato et 
al. 2000). 

For a given star formation history, gas fraction and 
metallicity, one of the parameters that largely affects the 
FIR emission is the escape time t esc . Silva et al. (1998) 
and Granato et al. (2000) have shown that the UV, opti- 
cal and FIR properties of local spirals are well reproduced 
with a typical t esc ~ 3 Myr. Normal star forming regions 
are also characterized by a moderate visual attenuation 
Ay ~ 1 mag. In the case of dusty starbursts, the slope 
of the UV continuum (Meurer et al. 1999) indicates large 
obscuration times. Silva et al. (1998) and Granato et al. 
(2000) were able to reproduce the SEDs, from the UV to 
the FIR, with a t esc in excess of several Myr and a com- 
pact geometry with a characteristic radius of a fraction 
of a Kpc. Similar obscuration times seem to be required 
to interpret the optical spectra of very luminous infrared 
galaxies (Poggianti et al. 2001). 

We have considered as representative of a normal spiral 
a model with the parameters shown in Table |l]. We adopt 
a Schmidt type star formation law, i.e. SFR=i/Mg, with 
infall of gas with primordial composition in a time scale 
ti n f, and a Salpeter initial mass function (IMF) (slope 
x=1.35 in mass), from 0.15 M Q to 120 M Q . For starburst 
galaxies we superimpose an exponentially decreasing burst 
(ib=10, 15, 25, 50 Myr) of star formation at an epoch of 
11.95 Gyr, with the same IMF. The mass of stars formed 
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Table 1. Parameters for a disk galaxy model. See text for 
details. 
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during the burst is set to 10% of the disk mass. Our main 
conclusions are not affected by details in these choices. 



3. Radio emission 

A clear picture of radio emission from normal galaxies is 
still missing, particularly for the often dominant non ther- 
mal component. Indeed, it is well known that the intensity 
of the thermal component is tightly related to the number 



of H ionizing photons, Q(H), and scales as i 



(Rubin, 



1968). Computing HII region models with CLOUDY (Fer- 
land 1996), for different mass (1O 4 M to 10 5 M Q ), metal- 
licity (Z = 0.008, 0.02 and 0.05) and age (1 to 10 Myr) 
of the ionizing cluster, and for electron densities of 10 1 to 
10 3 cm -3 and inner radii of the nebula of 10 and 100 pc, 
we have obtained the following average relation at 1.49 
GHz (see Panuzzo et al. 2002 for details of the inclusion 
of nebular emission in GRASIL): 



Q(H) ( T e \ 



0.45 



, — ( v Y d) 

erg/s/Hz 5.495 x 10 25 V 10 4 K/ V1.49GHz/ w 
It is then straightforward to obtain Q(H) and the intensity 
of thermal radio emission from simple stellar populations. 

On the other hand, little is known about the source of 
the non thermal emission which, in normal star forming 
galaxies may account 90% of the radio emission (Condon 
1992). Observations indicate that FIR and radio emission 
are strongly correlated over a wide range of IR luminosi- 
ties. At 1.49 GHz (Sanders & Mirabel 1996): 
F FIR /(3.75 x 10 12 Hz) 



q = log ■ 



2.35 ±0.2 



(2) 



i^(1.49GHz)/(W m^Hz" 1 ) 
where F FIR = 1.26 10- 14 (2.58 5 60/1 m + Sioo„m) W m" 2 , 
with Sgo and S'ioo in Jy This correlation suggests that 
NT emission is related to the recent star formation and the 
most likely mechanism is synchrotron emission from rela- 
tivistic electrons accelerated into the shocked interstellar 
medium, following CCSN explosions. But the poor knowl- 
edge of the accelerating mechanism hinders any quantita- 
tive prediction of this phenomenon which can only be for- 
mulated in an empirical way (Condon 1992). The observed 
average luminosity per supernova event can be estimated 
with the ratio between NT radio emission and CCSN rate 
{vccsn) in our Galaxy (e.g. Condon & Yin 1990). With 
vccsn — 0.015 (Turatto, private communication) and 
after converting L Q 4 qjj z — 6.1 x 10 21 W Hz^ 1 (Berkhui- 
jsen, 1984) to 1.49 GHz by assuming a radio slope of 
a=-i" 



- d }° gS " =0.8 1 we obtain 

a log v 



vccsn/y? 1 



(3) 



It is still a matter of debate why the FIR/radio corre- 
lation is unaffected by the dependence of the electrons life- 
time and luminosity on the magnetic and radiation density 
fields, which may change significantly in different environ- 
ments. 

Lisenfeld et al. (1996) claimed that, in normal star 
forming galaxies, the existence of the FIR/Radio corre- 
lation requires a correlation between the radiation and 
magnetic energy density fields. However the radiation field 
changes dramatically in starburst galaxies, and there must 
be a significant fine tuning between the parameters reg- 
ulating the intensity of the two fields, because they must 
scale in such a way that they give rise to the FIR/Radio 
correlation and, at the same time, they must prevent in- 
verse Compton to dominate over synchrotron losses. 

As a possible way out of this conundrum, we suggest 
here that a FIR/Radio correlation originates because syn- 
chrotron electron lifetimes are shorter than the fading time 
of the CCSN rate. Assuming that cosmic ray electrons 
are injected during the adiabatic phase of SN explosions 
(t<10 4 -10 5 yr), i.e. in a characteristic time scale which is 
short compared to the star formation time scale, the bolo- 
metric synchrotron luminosity at an epoch T is given by 

jNT _ / ' VccSN{T -t)l NT {t)dt (4) 
Jo 

where t is the lookback time, l NT (t) is the NT luminosity 
of the injected electrons after a time t has elapsed, T is 
the age of the galaxy and rf is the lifetime of electrons 
against synchrotron losses and depends on the intensity 
of the magnetic field. For normal spiral galaxies rf (<< 
T) may be as large as several 10 7 yr but, since the star 
formation rate is almost constant over the last Gyr, the 
above integral becomes 



- NT 



vccsn 



o 

where we put l NT (t) = dE/dt. Thus the bolometric NT 
luminosity scales linearly with the SN rate and the propor- 
tionality constant is the injected energy of the electrons 
per SN (E ei ). 

In the case of starburst galaxies the SN rate cannot 
be considered constant over such a large time scale. How- 
ever also rf must be much shorter: to avoid significant 
losses from inverse Compton on the intense stellar radia- 
tion field rf « 1-0.1 Myr (Condon 1992). Thus we may 
still make use of the approximation in Eq. |H[ because we 
may consider vccsn almost constant over such a small 
time scale. 

In brief, the NT radio luminosity of a galaxy is pro- 
portional to the integral of the synchrotron power over 
the electron lifetime, and an increase of the former in a 
larger magnetic field is compensated by a shortening of 
the latter. 

Since both the SN rate and the FIR emission are 
strictly related to the recent star formation rate, our jus- 
tification of the validity of Eq. ^| in very different envi- 



dE , 
— dt 
dt 



vccsNE e 
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ronments, explains why the FIR/Radio correlation is so 
robust. No fine tuning is necessary, apart from the require- 
ment that the magnetic field is large in star burst galaxies, 
which, by itself, is an independent observational fact. 

The considerations above apply to the bolometric radio 
luminosity, but what is actually measured is the specific 
luminosity L„ T . However, under plausible assumptions, 
also L^ T scales linearly with the SN rate, with the only 
dependence on environmental conditions being an almost 
vanishing one on the magnetic field. This is a consequence 
of the observed spectral index of the NT radio emission. 

Indeed, the diffusion-loss equation (e.g. Longair 1994) 
for the time evolution of the number density of relativistic 
electrons per unit energy interval N(E,t) is 

= ±[ h {E)N{E,t)]+Q{E,t)+DW 2 N{E,t){&) 

where b(E) is the rate at which particles lose energy, 
Q(E,t) is the rate at which electrons are injected in the 
system per unit volume, time and energy interval (the 
source term) and D is the diffusion coefficient. We assume 
that the distribution of sources of fresh electrons (i.e. the 
SNae) is sufficiently uniform and extended to make neg- 
ligible the diffusion term, and that each SN event injects 
high energy electrons with a power law energy spectrum 
cx E~ p , so that the source term can be written as 

Q{E,t) = k UCCS v N{t) E-P 

where k and p are constants (depending on the detailed 
physics of SN explosion) and V is the volume of the 
system. Also, we are interested in stationary solutions 
dN/dt — 0, that is we consider electron energies at which 
the lifetime against radiative losses is short compared to 
the typical time scale for variations of vccsn^)- The so- 
lution of Eq. H is then simply 
N(F] kv CCSN E-^) 
N{E) ~ V( P -l)b(E) 

where b(E) is the rate at which particles lose energy. If 
the dominant loss mechanism is synchrotron radiation, 



b(E) 



dE 



,B Z 



— = -<r T cy — 

at synch O 07T 



In this case 

N(E)= V . B 2 

where the multiplicative factor A depends only on p and 
fundamental constants: electrons with a sufficiently short 
lifetime have a power law energy distribution, with an 
index steeper by 1 than the injected distribution. 

On the other hand, the synchrotron luminosity of an 
optically thin source with a random magnetic field and 
electrons having a power law energy distribution N(E) — 
CE-i is 

L^ T = ju,synchV = CD{q) V 

D(q) depends only on q and fundamental constants. This 
equation with C = AkvccsN /V B 2 and a — (q — l)/2 
becomes 

L? T = vccsn (AkD(p) B a -^- a ) 



which shows that also the specific luminosity is propor- 
tional to the SN rate, with a proportionality factor which, 
given the observed spectral index of radio emission a w 1 , 
has only a very weak dependence on the magnetic field. 
It may be useful to note explicitly that the above formula 
integrated over v yields again the result expressed by Eq. 
[|, namely that the bolometric NT luminosity does not 
depend on B. This is because the limit of integration in 
frequency are those corresponding to the limits in energy 
within which the power law electron distribution applies, 
and the link between the two is v cx B x E. 

The crucial assumptions in this derivation are that the 
SN rate is constant over time scales of the order of the life- 
times of relevant electrons, and that the dominant electron 
energy loss is synchrotron emission, as already discussed. 



3.1. The contribution of Radio Supernova Remnants 

Before discussing how to calibrate Eq. |[ we consider the 
possible sources of non thermal radio emission in our own 
Galaxy. 

We have checked that, among the identified galactic 
sources, only radio Supernova Remnants (SNR) may pro- 
vide a significant contribution to the NT radio emission. 
Other sources like the pulsars themselves and the bubbles 
of NT radio emission associated with X-ray binaries pro- 
vide a negligible contribution. However SNRs, while being 
the obvious most appealing sources, cannot be responsible 
of the bulk of NT radio emission in normal galaxies for two 
reasons. First, their spectrum has a characteristic radio 
slope L v cx iy a with a broad range of a between 0.2 and 
0.5 (Gordon et al. 1999) which is, on the average, less than 
that characteristic of normal galaxies (a ~0.8). Second, 
one may easily show that they cannot supply more than 
about 5% of the total NT luminosity, or, equivalently, the 
SN rate needed to reproduce the NT radio emission of the 
Galaxy is more than one order of magnitude larger than 
observed (see e.g. Condon 1992 and references therein). 

In fact, because the typical lifetime of a SNR is of a 
few 10 4 yr, we may estimate the contribution of a popula- 
tion of SNRs originated by an instantaneous burst of star 
formation, by applying the fuel consumption theorem of 
post main sequence evolutionary phases (Renzini & Buz- 
zoni 1986). Indeed short evolutionary phases beyond the 
main sequence provide an integrated luminosity 

dm d '" 



L 



SSP 



dt 



ld(t)dt 



(7) 



where md is the dying mass (~ Mrum-off), ld(t) is the 
luminosity evolution within the phase, (j>(md) is the IMF, 
is the time derivative of the dying mass, and the inte- 
gral refers to the particular phase. Altogether the quantity 
4>{ m d)^ d ^ L provides the evolutionary flux, namely the rate 
of dying stars of the given population. 
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The luminosity evolution of a single SNR may be 
derived by combining the observed surface brightness- 
diameter (E - D) relation at 408 MHz 
E (W m _2 HsB _1 sr- 1 ) ~ 10" 15 (8) 
with the time evolution of the linear diameter of a SNR 
(Clark & Caswell 1976) 

D^-OASE^n- 1 / 5 ^ 5 (9) 
It follows that 

/(erg s^Hz- 1 ) ~ 3.55 x 10 26 D~* (10) 
and the evolution of the luminosity of a SNR (at 1.49 GHz, 
assuming a slope a=0.3) is (see e.g. Condon 1992) 

erg/s/Hz 

The integrated contribution of the population of SNRs 
from an instantaneous burst of star formation, at an age 
t, is then 

<Km d )^ (12) 



(11) 



tSSP 
-^1.49 



1.7 x lO 26 ^ 



-1/5 l/5_3/5, 



erg/s/Hz " ~" ~ ou ' c " rv "" v dt 

where E$q is the SN blast energy in units of 10 50 erg, 
n (cm ) is the ambient particle density and r c is the 
lifetime of the SNR in yr. Note that is a function 
of time, but it is almost constant during the considered 
phase. 

The lifetime of a SNR is usually associated with its 
adiabatic phase (Condon 1992) 



2 x 10 4 Ex 



o() 



yr 



thus 



SSP 
49 



0.6xl0 29 n- 2 /% 1/17 ^W 



dm a 



(13) 



(14) 



erg/s/Hz &u rv ' dt 

Integrating over the past SFR, ip(t), until we still get 
CCSN events (t < tccSN) and after defining the current 
supernova rate 

r-tccsN drrid 
vccsN(y*~ ) = / <f>(m d )——<p(t)dt. (15) 



jo dt 
we obtain for the average SNR non thermal luminosity per 
supernova event 
SNR Lf£/7(10 30 erg/s/Hz) 



0.06 n^ 17 E- 1/17 (16) 



77 

& 1.49 — / _1 

vccs N/yr 

This result, independent from the assumed IMF and star 
formation history, is also not very sensitive to the envi- 
ronment. It clearly shows that the contribution of SNRs 
is much less than required by Eq. || for plausible values of 
n and E$o (Condon 1992). In principle the lifetime of a 
SNR is longer than the adiabatic time and could also de- 
pend on the ambient density. Kafatos et al. (1980) quote 
for a hot cavity SNR tsnr = 1.4 x 10 5 yr, at which time 
the linear diameter is about 100 pc and the velocity of 
the shock front is 300 km/s, still capable of maintaining a 
shock temperature of about 10 6 K. However this lifetime is 
still one order of magnitude less than that required by the 
observed NT luminosity (L oc t 3 / 5 ) and, more important, 
when the diameter is larger than about 20 pc the observed 
E — D relation steepens considerably (E oc D~ 10 ) so that, 
in any case, the contribution beyond t c is negligible. 



4. Calibration of the NT radio emission 

In summary in modelling the properties of the NT radio 
luminosity in our Galaxy we have found that the only non 
negligible discrete sources SNRs can provide at most 6% 
of the NT radio luminosity (e.g. Condon 1992). Following 
Condon & Yin (1990) we have thus calibrated Eq. || (after 
accounting for the small contribution of SNRs) against 
the SN rate and synchrotron luminosity of our Galaxy. 
We assume for the NT radio emission 
L NT {v) 



10 30 erg s-iHz" 1 

/ zpSNRf V _ \-0.5 , jpel I V _ \-a\ v 

(£1.49(7^9) +Sl - 49( L49 ) ) X 



(17) 



VCCSN 

yr -i 



Since the NT contribution of SNRs has been evaluated in 
Eq. H {Ef% R ~ 0.06) then Eq. | requires that Ef A9 ~ 
1.38. In order to reproduce the observed average slope of 
the NT radio emission in normal spiral galaxies (~0.8), 
we have adopted a ~0.9. 

An independent check of this calibration is provided 
by the observed ratio between the FIR and radio emission 
in normal spirals (Eq. ||) . We have evolved several models 
with GRASIL for different values of the critical param- 
eters for the chemical evolution (different age, gas infall 
time scale and star formation efficiency), the escape time 
from molecular clouds (2 Myr and 3 Myr, Granato et al. 
2000) and the scale length of the dust distribution (3 Kpc 
and 4 Kpc ). The values of q turn out to be quite inde- 
pendent from the adopted parameters and cluster around 
q=2.35, in excellent agreement with the observations. We 
emphasize that, in our models, the SFR depends linearly 
on the gas fraction (Schmidt law), while the FIR emission 
depends on the SFR, gas fraction and metallicity, and the 
SN rate on the recent SFR. The consistency between FIR 
emission, radio emission and supernova rate, is thus re- 
markable and should be considered as a successful test of 
the model. We have also obtained the following calibra- 
tions between radio emission and star formation rate: 
SFi?(M Q /yr) = 2.8 10" 28 L(8.4GHz)/(erg/s/Hz) 
SFi?(M Q /yr) = 7.5 10~ 29 L(1.49GHz)/(erg/s/Hz) 
These relations hold only for the case of a constant SFR 
for more than 100 Myr and are in good agreement with 
those obtained by Carilli (2000). It is worth noticing that 
with a SN rate of 0.023 - the value adopted by Condon 
1992 - we would have obtained a SFR calibration higher 
by about 50%. 

Before concluding this section it is worth commenting 
on our choice of the lower limit of the initial mass giv- 
ing rise to type II SNe. It is commonly assumed that it 
corresponds to an initial mass of about 8 M Q for stan- 
dard convection while it lowers to about 6 M Q for models 
with convective overshoot. However, recent investigations 
on the evolution between 8 and 10 M Q (without overshoot) 
indicate that the final fate of stars born in this mass range 
is that of a white dwarf, instead of an electron capture SN 
(Ritossa, Garcia-Berro & Iben 1996). The same fate is ex- 
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Table 2. Starburst parameters for the SEDs in Figg. |] 
and | 



Case 


Age 


t b 


Current SFR 


Burst Mass 




(Myr) 


(Myr) 


(M /yr) 


(M/MtoO 


M82 










a dotted 


45 


50 


3.41 


0.011 


b solid 


24 


8 


2.37 


0.011 


ARP220 










a solid 


50 


50 


271 


0.121 


b dotted 


40 


30 


280 


0.121 


c dashed 


25 


8 


149 


0.121 




Fig. 1. M82: two models with different age and e- folding 
time scale (Tab. ||), fit the observed SED from UV to the 
sub-mm. Inclusion of radio wavelengths is necessary to 
disentangle the two cases. For details on the models see 
Silva et al. (1998). 

pected in the mass range 6 to 8 M Q if the overshoot scheme 
is adopted (Portinari, Chiosi, & Bressan 1998). We have 
thus assumed that CCSN are produced in stars with mass 
M>8 M© and for ages younger than tccSN — 50Myr. 

5. FIR/Radio properties of starburst models 

In this section we examine the FIR and radio properties 
of our new models and compare them with two well stud- 
ied local starbursts, M82 and ARP220. To obtain a realis- 
tic starburst model, we added an exponentially decreasing 
burst of star formation to the underlying secular disk evo- 
lution of the spirals models used in the previous section. 
Table ^| summarizes the burst parameters in few adopted 
models. Notice that in all models of the same galaxy, the 
total mass of gas converted into stars is the same (column 
M/Mtot)- The corresponding GRASIL parameters have 
been selected according to Silva et al. (1998) and, in par- 
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Fig. 2. Same as in Fig. |l|, for three models for ARP220. 

ticular, the adopted obscuration time must be significantly 
larger than that characteristic of normal galaxies. 

Fig. |l] shows two fits to the SED of M82, differing only 
in the age (45 and 24 Myr) and e-folding time (50 and 8 
Myr) of the burst. The UV-optical-FIR SED of the star- 
burst is degenerate, in the sense that it may be fitted by 
different set of burst parameters. Both models provide a 
good fit to the UV to sub-mm data. At radio frequencies 
however, only the short burst is able to reproduce the ob- 
servations while the longer one falls significant below the 
observed flux. A similar example is shown in Fig. |^ for the 
ULIRG ARP220. ARP220 is not consistent with a short 
burst, though this could not be excluded by the inspec- 
tion of UV-optical-FIR data alone. Notice that the current 
SFR of the models differ by about 40% and 90% in M82 
and ARP220, respectively. 

The cases of M82 and ARP220 show that the inclusion 
of the radio emission constitutes a powerful diagnostic tool 
to investigate starburst galaxies. Indeed, while the UV, the 
FIR and the radio thermal continua are sensitive to the 
number of living massive stars, the NT emission is a mea- 
sure of the current CCSN rate. Thus the FIR/Radio ratio 
is a measure of the ratio between the almost instantaneous 
SFR and the SFR averaged over the last few tens of Myrs. 
Combining the FIR and radio spectral regions is particu- 
larly important for the case of obscured starbursts, where 
the burst properties cannot be derived by UV, optical, 
NIR and even mid and far IR (continuum) observations. 
We ought to stress in the following, that a careful combina- 
tion of FIR and Radio data may be sufficient to determine 
the recent star formation history even for those galaxies 
that are not as thoroughly observed as M82 or ARP220. 
Inspection of Figg. |l| and || shows that the expected value 
of the q ratio at 1.49 GHz and 8.44 GHz changes by more 
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Fig. 3. The run of the star formation rate in our starburst 
models for different e-folding time scales, t& = 10, 15, 25 
and 50 Myr. 
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Fig. 4. Evolution of the ratio between the instantaneous 
star formation rate (SFR, M©/yr) and FIR luminosity 
(erg/s), see text for the definition. 



than a factor of 2 for the different models. This suggests 
that deviations from the average FIR/Radio correlations 
should be expected among obscured starbursts and actu- 
ally could be used as a powerful diagnostic for the analysis 
of the star formation in the burst. Furthermore, also the 
radio slope is affected because at younger ages the rel- 
ative contribution of thermal emission is larger than at 
older ages. 

To highlight this point we have carefully analysed the 
evolution of selected starburst models that possibly en- 
compass different realistic scenarios. For the star forma- 
tion rate during the burst we have assumed four different 
e-folding time scales, 4=10, 15, 25 and 50 Myr. The mass 
of stars formed during the burst was set to 10% of the un- 
derlying disk mass. The obscuration time t esc was set to 
linearly decrease with time, from 4 down to a minimum 
of 3 Myr, characteristic of normal galaxies. The run of the 
SFR of the models during the burst is depicted in Fig. |[ 

Fig. [| shows the ratio between the current SFR 
(M Q /yr) and the FIR luminosity, L FIR . The latter, repre- 
senting a fair measure of the luminosity between 40 and 
120 /mi, is defined in analogy with Helou et al. (1988) as 
L Fm (erg/s)= 1.257 x 10" 11 (2.58 L 60 + L wo ) where L 60 
and Lioo are the luminosities in erg/s/Hz at 60 and 100 
/itm. In this figure, the model with 4=25 Myr has been re- 
computed with a fixed t esc =3 Myr, corresponding to that 
of normal galaxies. Fig. |B| shows the ratio between the SFR 
(M Q /yr) and the radio luminosity (erg/s/Hz) at 1.49 GHz 
(upper panel) and at 8.44 GHz (lower panel). Obviously, 
in the latter figure, t esc is not relevant. Three main evo- 
lutionary phases can be recognised, the starburst, post- 
starburst and quiescent phase, characterized by different 
relations between either the FIR or the radio emission 
and the current SFR. At early times, during the starburst 
phase, the SFR/Flux ratio is significantly higher than the 
asymptotic value defined by normal star forming galax- 



ies. At later times, when the ratio falls below the value 
of normal galaxies, the models enter the "post-starburst" 
phase. Finally, as the burst extinguishes, the ratio turns 
to the asymptotic value defined for normal star forming 
galaxies. Depending on the evolutionary status and the e- 
folding time of the SFR, the relations between infrared or 
radio emission and SFR may change by even one order of 
magnitude. As already anticipated, NT radio emission is 
more sensitive to the past star formation history than the 
FIR emission and the corresponding ratio shows a larger 
variation. 

6. FIR/Radio properties of compact ULIRGs 

We have already anticipated that our models predict the 
existence of a correlation between the radio slope and 
the FIR/radio ratio, because both quantities are affected 
by the relative contributions of FIR, thermal and non- 
thermal radio emission that change during the evolution 
of the burst. Understanding whether the use of such cor- 
relation can provide further physical insight on the star 
formation process becomes, obviously, particularly rele- 
vant in the case of highly obscured starbursts such as 
ULIRGs, emitting the bulk of their radiation in the mid 
and far infrared. ULIRGs are characterized by extreme 
FIR luminosities exceeding the UV- optical power by even 
more than one order of magnitude, and are believed to 
be transient phases of galaxy activity associated with the 
dynamical interaction and merging of gas rich systems. 
Whether their powerful emission is of starburst origin or 
is due to the AGN is still under debate (Sanders et al. 
1988, Condon et al 1991, Veilleux et al 1994) but, re- 
cently, it has become clear that these objects may har- 
bour huge compact star forming complexes and the AGN, 
at the same time (Rowan- Robinson, 1995). However, it re- 
mains still unclear what fraction of the bolometric lumi- 
nosity is eventually provided by the central monster and at 
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Table 3. Selected galaxies from the sample of compact ULIRGs by Condon et al. (1991) 



Name 


Sl.49Gflz 


q 


S8.44Gflz 


Ql. 49-8. 44 


FIR 




Symbol" 


Classification* 


IRAS 10566+2448 


46.1 


2.55 


14.1 


0.68 


11.90 


- 


H 


HII 


A 1 1 m n_i_/i i ri7 


zo.u 


Z.04 


in 7 

1U. / 


U.oo 


11. OZ 




TI 
il 


TTTT 
Jill 


IRAS12112+0305 


22.6 


2.66 


10.0 


0.47 


12.18 


_ 


H 


HII 


UGC08335 


51.2 


2.46 


17.0 


0.64 


11.62 


590 


H 


HII 


IRAS22491-1808 


6.1 


3.00 


3.0 


0.41 


12.02 


- 


H 


HII Tb=5.2 


IRAS17132+5313 


28.4 


2.46 


8.9 b 


0.67 


11.79 


450 


H 


HII Tb<5 


UGC04881 


29.0 


2.50 


8.8 


0.69 


11.61 


210 


H 


HII Tb>7 a 


IRAS01173+1405 


43.1 


2.46 


12.7 


0.70 


11.54 


160 


H 


HII Tb>7 a 


Mrk331 


67.5 


2.51 


21.5 


0.66 


11.27 


510 


H 


TTTT rni r-r 

Hll 1 b>7 a 


IRAS10173+0828 


8.8 


2.92 


5.0 


0.28 


11.70 




H 




ARP220 


301.1 


2.63 


148.0 


0.41 


12.11 


- 


P 


HII Tb>7 a 


IRAS04191-1855 


27.3 


2.49 


8.9 


0.65 


11.34 


100 


L 


Liner Tb<5 


IRAS08572+3915 


6.5 


3.11 


4.1 


0.27 


11.96 


- 


L 


Liner Tb<5 


IRAS14348-1447 


33.2 


2.38 


9.7 


0.71 


12.17 


- 


L 


Liner Tb<5 


IRAS01364-1042 


17.0 


2.67 


8.2 


0.42 


11.67 


- 


L 


Liner Tb<5 


TTrrinoo'7 

UGC08387 


106.0 


2.28 


34.9 


0.64 


11.51 


1220 


L 


Liner 


Mrk273 


130.0 


2.31 


43.5 


0.63 


12.04 




L 


Liner Tb>7 a 


IRAS03359+1523 


18.9 


2.58 


11.0 


0.31 


11.37 


370 


A 


HII Tb>7 


Mrk848 


46.8 


2.38 


12.1 b 


0.78 


11.72 


240 


A 


HII Tb>7 


UGC02369 


42.7 


2.39 


13.3 


0.67 


11.42 




A 


HII Tb>7 1 


IIIZw035 


39.3 


2.58 


19.7 


0.40 


11.46 


590 


A 


Liner Tb=5.4 


IRAS15250+3608 


12.8 


2.81 


10.5 


0.11 


11.88 




A 


Liner Tb=6 


UGCObright 


146.0 


2.10 


52.6 


0.59 


11.93 




A 


Liner Tb>7 1 


NGC2623 


97.8 


2.51 


35.5 


0.58 


11.47 




A 


AGN Tb>7 1 


NGC0034 


58.7 


2.52 


15.2 


0.78 


11.28 




S 


Sey2 Tb=5 


Zw475.056 


26.0 


2.65 


8.2 


0.67 


11.37 


190 


S 


Sey2 Tb<5 


IRAS05189-2524 


28.1 


2.76 


11.4 


0.52 


11.91 




s 


Sey2 Tb<5 


Mrk231 


240.0 


2.24 


265.0 


-0.06 


12.35 




M 


Seyl Tb>7 



b Uncertain value 



* Electron density, cm -3 , from Veilleux et al. (1995) 
^ Symbols adopted in Figg |^ and ^ 

* HII, Liner, AGN, Seyfert based on Veilleux et al. (1995); TJ, logarithm of the brightness temperature from VLBI data by 
Smith, Lonsdale & Lonsdale (1998) 



what level this affects FIR and radio emission. For exam- 
ple, very recently Berta et al. (2002) and have shown that 
in IRAS 19254-7245 ("The Superantennae"), the AGN is 
contributing about 40-50% of the MIR and FIR emission 
( essentially the bolometric luminosity), while in the ra- 
dio it overcomes by one order of magnitude the emission 
from the starburst. At the same time, for IRAS 20100- 
4156, one of the brightest nearby objects, both the multi- 
wavelegth analysis (Fritz et al 2002), spectro-polarimetry 
(Pcrncchele et al 2002) and near-IR spectroscopy (FWHM 
Paa ~440Km s _1 , Valdes et al. 2002) show no evidence 
of the AGN. 

In this respect, the sample of compact ULIRGs se- 
lected by Condon et al. (1991) turns out to be particu- 
larly interesting because high resolution VLA maps sug- 
gest that the bulk of their radio continuum could be of 
starburst origin and several objects have been observed 



at both 1.49GHz and 8.44GHz, so that a direct compari- 
son in term of " observables" can be made with our star- 
burst models. We thus isolated all the compact ULIRGs of 
the Condon et al. sample observed at both 1.49GHz and 
8.44GHz, in Table ||. Optical spectroscopic classification 
by Veilleux et al (1995) indeed showed that a number of 
the selected objects have emission line ratio characteris- 
tic of photoionization by massive stars (HII galaxies) , but 
other show higher excitation and were classified as Liners 
or even Seyfert 2 galaxies. Mrk 231 is a Seyfert 1 object. A 
subset of the original Condon et al. sample has been sub- 
sequently mapped with VLBI by Smith, Lonsdale & Lons- 
dale (1998, SLL), to establish the nature of the milli-arc 
seconds structures and in particular to obtain the bright- 
ness temperature of the compact inner core. SLL did not 
found significant correlations between the VLBI emission 
(usually 10% of the 1.6GHz total flux density) and other 
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Fig. 5. Evolution of the ratio between the instantaneous 
star formation rate (SFR, Mq/jt) and radio luminosity 
(erg/s/Hz) at 1.49 GHz (upper panel) and at 8.44 GHz 
(lower panel) 



physical parameters like total radio power, FIR emission 
and radio slope. Perhaps more important, they have found 
that that optical excitation is not strictly correlated with 
high T5. As can be seen from their data, summarized to- 
gether with the Veilleux ct al. classification, in the last 
column of Table 0, there are HII objects showing high 
cores and LINERs and Sey 2 (all) galaxies with central 
T5 consistent with a starburst origin. Furthermore, SLL 
were able to show that some of the high Tj, cores could 
be explained by bright radio supernovac complexes, con- 
sistent with the star formation rates required by the FIR. 
These objects are marked with a lower case "a" in the 
last column of Table |L while those that could not be in- 
terpreted in terms of bright SNRs by SLL, are marked 
by a lower case "1". Some of the high T& objects had not 
enough detailed structure to allow this analysis. 

According to the information provided by Veilleux et 
al. (1995) and SLL, we grouped all objects into four broad 
categories, named with the symbol reported in column 8 
of Table |. The letter "H" indicates all galaxies of HII 
type, but for three objects for which SLL could not estab- 



lish that their high TJ, is compatible with RSN complexes. 
Analogously, Liners are indicated by the symbol " L" when 
their T& is missing or is low, or their high Tb was found 
compatible with RSN complexes. Object classified as HII, 
Liners or AGN that do not fulfil the above requirements 
are indicated by the symbol " A" while, for Seyfert 2 galax- 
ies we use the symbol "S". Finally "P" and "M" refer to 
ARP220, a typical starburst galaxy and Mrk 231, a Seyfert 
1 galaxy, respectively. According to this grouping, the ef- 
fect of the AGN, if any, should not be significant in types 
H to L. For example a recent comprehensive analysis of 
optical-NIR-FIR SED and of spectro-polarimetric data of 
IRAS 22491-1808, indicates that if the AGN is present in 
this object, it must contribute less than 5% of the total 
FIR emission (Fritz et al 2002). 

In Figure |^ we show the run of the q values expected 
from our starburst models, Eq.(||), at 1.49 and 8.44 GHz, 
against the FIR luminosity (now in solar units). Super- 
imposed are the selected compact ULIRGs from the local 
sample of Condon et al. (1991), with the corresponding 
symbols listed in Table |[ 

The models can be arbitrarily shifted in the horizon- 
tal direction by changing the mass of stars formed in the 
burst which, for the adopted normalization, corresponds 
to about 10 10 Mq. To fit the IR luminosity of the brightest 
compact sources we need a SFR higher by about a factor 
of 2 to 3, for the shorter burst. We note that the observed 
dispersion of the q ratio is well explained by the starburst 
models, in particular in the lower panel, where, as we will 
see, the effects of free-free absorption should be negligible. 
The mass of the burst is however not known a priori and 
while Fig. || shows that the observed data are consistent 
with the sources being compact starbursts, it does not al- 
low to obtain an estimate of the burst parameters. On the 
other hand, the case of M82 and ARP220 (Figg. [I] and [§) 
indicates that this is within the possibility of the model 
when enough observations are available. 

We thus combine in Fig. |?| the FIR/Radio ratio with 
the slope of the radio emission a between 1.49 and 8.4 Ghz. 
Both quantities are independent from the intensity of the 
burst. Rather, they depend on the form of the recent star 
formation history, so that the path in this diagram traces 
the evolutionary status of the starburst. The radio slope of 
the models changes because of the variation of the domi- 
nant source of radio emission as the starburst ages. During 
the first 3 to 4 Myrs only thermal emission from HII re- 
gions contributes to the radio emission, and the radio flux 
has a characteristic slope a ~0.1. Then CCSN explosions 
feed relativistic electrons into the galactic magnetic field, 
and NT emission steepens the spectrum toward a slope 
which is more typical of normal galaxies. At the same time 
the total radio power increases and the evolution of the q 
ratio, though affected by the corresponding increase of the 
FIR emission, continuously decrease to a minimum value. 
At this stage the SFR has decreased significantly and the 
model can be considered in a post starburst phase. At 
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Fig. 6. Evolution of starburst models in the q vs Lfir 
diagram. The upper panel is for q at 1.49 GHz, the lower 
one is for q at 8.44 GHz. Stars indicate the positions of 
typical models of quiescent star forming disks. The model 
with £(,=25 Myr and t esc =3 Myr is also shown (small solid 
dots). Symbols refer to the subset of compact ULIRGs 
(Condon et al. 1991), drawn in Table [| 




a 



Fig. 7. Upper panel: evolution of the ratio q vs radio slope 
a between 8.4 and 1.4 GHz for our starburst models. Char- 
acteristic starburst ages range from log t (yr) = 6.3 to 8.6. 
The model with tb=25 Myr has been recomputed both 
with a t esc =3 Myr (crossing the observed data) and by 
adopting a free-free absorption with t v =i.&GHz = 1 (model 
on the left side up to its minimum value of q). Lower panel: 
same for z/=8.4 GHz. Note that q is no more affected by 
free-free absorption. Symbols refer to the subset of com- 
pact ULIRGs (Condon ct al. 1991), drawn in Table §. 



even older ages the ratio increases again, but this corre- 
sponds to the very late phase where the model can no 
longer be considered representative of an ultra luminous 
galaxy. The model with i;,=25 Myr and obscuration time 
of normal galaxies, t esc =3 Myr, is not able to reproduce 
the observed high values of q, in both panels, confirming 
that the low UV flux observed in other obscured starburst 
(see e.g. Figg. [l] and ||) indicates that this escape time is 
too short. The observed data show a clear trend of increas- 
ing q at decreasing slope, in the 1.49 GHz plot, with the 
models being only able to delineate an upper envelope. 
At a higher frequency (lower panel), the trend in q disap- 
pears while the models still delineate the upper envelope. 
A possible explanation is that part of the effect is caused 
by free-free absorption (see also Condon et al. 1991). At 



radio frequencies the optical depth of free-free absorption 
by a cloud with electron density N e and size 1, is 

4* - 8.2 x 10- 2 T- 135 (-^)- 21 (^^) (18) 

{jtlz pc/cm 

In order to highlight the effects of free-free absorption we 

f f 

have recomputed a model with optical depth t( 49Ghz = l 
(dot-dashed line on the left side in Fig. |^). Assuming an 
electron temperature T e = 10 4 K, this optical depth cor- 
responds to an emission measure of about 6xl0 6 pc/cm 6 . 
The latter can be achieved by considering 2.5xl0 8 M Q 
within 250 pc or 7.8xl0 8 M Q within 390 pc of ionised gas, 
for an average electron density of about N e =150 cm -3 
or iV e =125 cm -3 and solar composition, respectively. The 
above figures compare fairly well with electron densities 
derived from emission lines ratio by Veilleux et al. (1995), 
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listed in Table ||. They are consistent both with the fact 
that the selected ULIRGs are compact radio sources with 
typical sizes less than 1 Kpc and possibly as low as a few 
hundred pc (Soifer et al. 2000), and that the total mass 
involved in our starburst models is of the order of 1O 1O M0. 

The distribution of the data in Fig. [?] also renders quite 
unlikely the possibility that the slope variation is due to 
an increasing importance of the electron cooling by in- 
verse Compton. Indeed it could be that in ULIRGs with 
the lower a electrons cool down by inverse Compton on the 
stellar radiation field: in those objects only thermal emis- 
sion would be present at radio wavelengths. Not only this 
would again require a fine tuning of the two cooling pro- 
cesses, in order to give rise to a distribution of slopes be- 
tween 0.1 and 1 but, and more important, it would be dif- 
ficult to explain the absence of the trend at the higher fre- 
quency. This enforces our interpretation in terms of free- 
free absorption because, given the frequency dependence 
(Eq. |l8|), this effect should not affect the value of q at 8.44 
GHz. The quantity q g 44 Qjjz' which still shows a range of 
about 0.6dexp, can be considered as a genuine measure of 
the age of the compact starbursts. Unfortunately, in this 
diagram, the slope is still affected by free-free absorption, 
and it is not possible from the above data alone to identify 
precisely the evolutionary status of the ULIRGs. However, 
our models suggest that, observing ultra luminous galax- 
ies at frequencies between 8.44 GHz and 23 GHz would 
possibly constitute a powerful tool to investigate on the 
recent star formation history of obscured starbursts. 

From the upper panel of Fig. [j] it is also evident that 
there is no sharp threshold value of ^ 49 qjj z that can 
be safely used to delineate a separation between star for- 
mation and AGN powered ULIRGs. There is a tendency 
for objects classified "A" to occupy the lower boundary 
allowed by starburst models but, even in the case of the 
Seyfert 1 galaxy Mrk 231, the FIR brightest object in 
our sample, it would be difficult to exclude a starburst 
origin resting only on its q 1 49 qh z ■ However its location 
in the corresponding q-slope diagram is not matched by 
any of our starburst models and, by looking at the higher 
frequency data q g 44 qjj z , it appears to be at least three 
times more radio powerful (relative to the FIR) than all 
the other sources and than allowed by our models in the 
starburst phase. Is this an indication that the parame- 
ters plotted in Fig. [t] (possibly making use of a higher 
frequency to avoid free-free contamination of the radio 
slope) could be used as a diagnostic diagram to disentan- 
gle a starburst from an AGN? And, in this case, why are 
all the Seyfert 2 galaxies in our sample just populating the 
upper envelope of starburst models?. From the above plot 
we may only conclude that a minimum value of q x 49 qjj z 
for a starburst powered source is around 1.8, and this may 
set a threshold value below which another source of radio 
emission has to be invoked. One should consider however 
that at this frequency the nuclear engine may be masked 
by free-free absorption. At 8.44 GHz free-free absorption 
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Fig. 8. The q ratio as a function of time for the models in 
Fig. @. 

should be less important and a fair threshold for star for- 
mation dominated ULIRGs is q g 44 GHz = ^- 

7. The post starburst phase 

In a recent comprehensive analysis of the FIR/Radio cor- 
relation in nearby Abell clusters, Miller & Owen (2001) 
found a statistically significant excess of star forming 
galaxies with enhanced radio emission relative to the FIR, 
toward the cluster centres. High resolution radio images 
have also excluded a significant AGN contribution to the 
radio emission in these galaxies. 

There is a long standing debate on the nature of this 
excess of galaxies with a low value of q 1 49 GHz > w hich 
is always interpreted as a radio enhancement and may 
reach a factor of three. Gavazzi & Jaffe (1986) advanced 
the hypothesis that the radio excess is caused by the ram 
pressure strengthening of the galaxy magnetic field, as the 
galaxy travels through the intracluster medium. Accord- 
ing to Miller & Owen (2001) this explanation seems not to 
work because of the lack of any correlation with the galaxy 
velocity. They suggest instead a compression of the galac- 
tic magnetic field by thermal pressure of the intracluster 
medium. 

However, due to its small scatter, the FIR/Radio cor- 
relation should be quite independent from the effects of 
the environment. We thus advance the hypothesis that 
the observed excess is due to an excess of post starburst 
galaxies in the central regions of the cluster. Indeed soon 
after a peaked starburst episode and/or a star formation 
interruption, the radio emission fades less rapidly than the 
FIR emission, causing an apparent radio enhancement. In 
this case the effect would be only indirectly due to the 
environment, as the gas rich galaxies enhance and/or ex- 
haust their star formation through a central crossing. 

All our models show a more or less pronounced min- 
imum value of q during the post starburst phase. The 
model with i;,=10 Myr remains below qi.49GHz=2 for 



A. Bressan et al.: FIR and Radio emission in dusty starbursts 



13 



about 20 Myr (Fig. g), while its luminosity at the min- 
imum q has decreased by one order of magnitude (Fig. 
but it is still infrared luminous. By decreasing the total 
mass formed in the burst it is easy to populate for a short 
time the region with qi. 49GHz below 2 and Lfir between 
10 10 and 10 11 L©. Thus a low value of qi.4GHz could simply 
be a natural consequence of the particular star formation 
history experienced by the galaxy. Because high resolu- 
tion radio images probe that the emission is extended, the 
recent star formation history, either simply interrupted or 
enhanced and exhausted in a burst, must have been glob- 
ally synchronised. Our models place an upper limit of less 
than 100 Myr to the age of the last major burst and/or 
interruption of the star formation. At a typical velocity of 
1000 Km/s, the galaxy has moved by only 100 Kpc since 
the beginning of the burst. This figure could be about two 
times larger if we allow for the formation of CCSN down to 
5 M Q (overshoot models without O, Ne, Mg white dwarfs) 
and use a slightly slower star formation decline. 

In summary, our models indicate that the star for- 
mation switch-off happened not too much far from the 
present galaxy position. The excess of low q galaxies is 
thus simply due to the larger probability of switching off 
the star formation in the higher density regions of the 
cluster. It is a nurture effect, possibly due to the higher 
degree of harassment suffered in the central regions. As a 
clear implication, our spectrophotometric models predict 
that these galaxies should show enhanced Balmer absorp- 
tion features, e.g. EW US >5A(see also Miller & Owen 
2001). If our interpretation is correct, the analysis of the 
statistics of the deviation from the FIR/Radio correlation 
at values of q lower than the average, coupled with the 
low characteristic time involved, provides an independent 
measure of the local rate at which the Butcher-Oemler 
effect is operating within galaxy clusters. 

8. High redshift starbursts 

In the last few years, a wealth of observations performed in 
FIR/sub-mm spectral regions revealed the existence of a 
new class of galaxies interpreted as the high- z analogue of 
the local ULIRGs (e.g. Smail, Ivison & Blain 1997, Hughes 
et al. 1998, Barger, Cowie, & Sanders 1999). The sub-mm 
fluxes, probably mostly powered by star formation rather 
than AGN, (Granato, Danese, & Franceschini, 1997; Al- 
maini, Lawrence, & Boyle, 1999), imply star formation 
rates of several hundreds of M /yr. 

The discovery of these galaxies has introduced a new 
test to the theories of structure formation, but to this aim 
the knowledge of their redshift distribution is of funda- 
mental importance. Due to the uncertain position, or to 
the lack, of the optical counterparts, a spectroscopic red- 
shift is available only for a few sub-mm sources. Instead, 
an estimate of the redshift for most of these galaxies has 
been performed exploiting the FIR-radio correlation ob- 
served for local star forming galaxies, under the hypothe- 



sis that high-z galaxies obey the same correlation as local 
ones (Carilli & Yun 1999, 2000; Dunne, Clements, & Eales 
2000; Yun & Carilli 2002). 

Due to the very different power laws of the sub- 
mm and the radio spectra, the spectral index sf 5 4 = 
log(5 , 3 53 /»S'i.4)/log(353/1.4), between the 850 /mi ' (353 
GHz) and the 1.4 GHz fluxes, depends very strongly on 
redshift, with the main uncertain factors in its effective- 
ness as a redshift indicator being (Carilli & Yun 2000; 
Dunne et al. 2000): (a) the dust temperature distribution, 
(b) the slope of the dust emissivity in the sub-mm and (c) 
the slope of the radio spectrum. Additional uncertainties 
are (d) the assumed validity of the FIR-radio correlation 
at high-z, that recently has been directly tested to z^ 1.3 
(Garrett 2002), (e) the possible presence of AGNs (whose 
values of sf 5 ^ could be attributed to lower z star forming 
galaxies rather than high-z AGN), and (f) the quench- 
ing of the synchrotron photons due to inverse Compton 
of relativistic electrons off the CMB radiation, which is 
expected to be important for z > 6 (Carilli & Yun 1999). 

With our model, we have investigated the (z) rela- 
tion for starburst galaxies, keeping into account the uncer- 
tainties a, b and c. These factors affect the intrinsic shape 
of the galaxy SED. Uncertainties due also to observational 
and calibration errors are accounted for by Hughes et al. 
(2001). The T distribution of dust and its slope are degen- 
erate in the sub-mm, because a shallower slope can mimic 
the effect of a colder dust component with a steeper slope 
(e.g. Silva 1999). This is due to the fact that the convo- 
lution of gray bodies of different temperatures (resulting 
from the T distribution of dust within a galaxy) yields a 
slope shallower than the effective one of the emissivity of 
dust. Moreover both the T distribution of dust in galax- 
ies and the slope of the radio emission depend on several 
factors that change during the evolution of the starburst 
(i.e. the evolution of the stars heating the dust and their 
distribution across the galaxy, the dust optical depth, the 
relative importance of free-free and synchrotron emission 
and the amount of free-free absorption, see Sect. |). 

Therefore we have adopted the set of starburst mod- 
els described in Sect. |[ i.e. the exponentially decreasing 
starbursts with e-folding times if, = 10, 15, 25, 50 Myr, 
complemented with one more case, tb = 5 Myr, and with 
the models for M82 and ARP220 (Figg. and |) , to study 
the dependence of sf 5 ^ (z) on the evolutionary status of 
the starbursts and the consequent different shapes of the 
SEDs. In particular, the model for ARP220 needs a slope 
for the emissivity of dust in the sub-mm of 1.5, while this 
is not required for M82 and other nearby galaxies (see 
Silva et al. 1998), whose SEDs can be reproduced with a 
slope of 2. In Sect. || we showed that the starburst models 
cover the range of observed values of the FIR-radio re- 
lation and of the radio spectral index. In Fig. [| we have 
checked that the range of values of the sub-mm to radio ra- 
tios (850/im— 1.4GHz and 450/im— 5GHz to check also the 
spectral regions observed at high z) of the models, cover 
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Fig. 9. Color-color plot log S S 50^ m / S 1A ghz vs 
log Sisopm/ Sschz- Squares: Starburst models with 
tf, = 5,10,15,25,50 Myr evolved to an age of 3^. Dia- 
monds: models for M82 and ARP220 (see Figg. @ and |), 
reproducing the SED from the UV to the radio (see Silva 
et al. 1998 for details). Asteriscs connected by solid line: 
evolution during the starburst of the models reproducing 
M82 and ARP220. Triangles: Data from Dunne et al. 
(2000), Condon et al. (1991). 



at least all the observed range. Furthermore, we can take 
into account also ratios not observed in the available local 
galaxy samples, but that might be expected to be present 
in some phases (e.g. during the first ~ 3 Myr since the 
onset of the starburst when the radio emission is mostly 
thermal). The models are evolved to an age of 3ib since 
the start of the burst. In Fig. [l(] we show the sf 5 4 vs z re- 
lation resulting by considering all our models (we plot the 
highest and the lowest one, corresponding respectively to 
a very initial phase dominated by free free thermal radio 
emission and a late phase, dominated by NT radio emis- 
sion), together with the relations by Carilli & Yun (2000) 
and Dunne et al. (2000). As already remarked by these 
authors, the effect of low frequency free-free absorption 
on the relation is important only at very low z, due to its 
strong frequency dependence (r/^ oc v~ 2A ). 

We note that the dispersion of local observations is en- 
tirely compatible with the one introduced by the evolution 
of the starburst. Furthermore, the different evolutionary 
stages and thus SED shapes of the models, result in a sig- 
nificant dispersion of s\ 5 f at any z, sufficient to invalidate 
any quantitative estimate of z. We show that it might be 
possible to reduce the uncertainty in the redshift estimate 
by combining sf 5 ^ with another spectral index that de- 
pends mainly on the starburst age. The latter is a radio 
spectral index (see Sect. ||and Fig. Indeed for redshift 
between and 6 the radio spectral index of each model 
changes less than 0.2. Thus one may use this index to con- 
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Fig. 10. The spectral index sf 5 4 as a function of redshit 
for our models (upper continuos and lower dotted lines), 
and from Carilli & Yun (2000) and Dunne et al. (2000) 
(respectively dashed and 3 dot-dashed lines for average 
and ±lcr). The upper model curve corresponds to a very 
initial phase of the starburst, when free free radio emis- 
sion dominates, while the lowest one is for a late phase. 
Data for galaxies and AGN with spectroscopic redshift 
are reported: Asteriscs are local galaxies by Carilli & Yun 
(2000) and Dunne et al. (2000); Diamonds are AGN by 
Rowan- Robinson et al. (1993), Ivison et al. (1998), Yun 
et al. (2000); Filled circles, Crosses, Squares, and Circles 
are sub-mm and radio selected galaxies from Barger et al. 
(2000), Smail et al. (2000), Eales et al. (2000), and Ivison 
et al. (2000) respectively. The Triangle is a cD galaxy from 
Edge et al. (1999). 

fine the sf^-z evolution within the possible values allowed 
by the corresponding selected starburst phase. 

In Fig. [ll] we show that a linear combination of the 
spectral index sf 5 4 with a radio spectral index (5 — 1.4 
GHz in the figures) may indeed reduce the redshift un- 
certainty at each z originating from age dispersion. As an 
example, Az at z = 2 for models is 3.12, 1.67, 0.82 respec- 
tively for Fig. [n] and the two panels of Fig. [ll|. On the 
other hand it is important to note that, by considering the 
few available sources (excluding the known AGN) with a 
spectroscopic z, the redshift range that would correspond 
to their ordinate value depends on the adopted spectral 
index or combination. A detailed investigation of the red- 
shift distribution of sub-mm galaxies is beyond the scope 
of this paper. 

9. Discussion and conclusions 

In this paper we have revisited the nature of the 
FIR/Radio correlation observed in star forming galaxies. 
To understand its origin and range of validity we have uti- 
lized models of normal star forming and starburst galax- 
ies. The infrared emission has been estimated with our 
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Fig. 11. Linear combinations of sf 5 % and a\ 4 for mod- 
els and data. For data with only one radio flux available, 
the radio spectral index was set to 0.7. See Fig. for 
references to data. 



population synthesis code, GRASIL, which is particularly 
suited for the prediction of the SED of star forming galax- 
ies, from the UV to the sub-mm. As for the radio emis- 
sion we have adopted a model which extends the capabili- 
ties of GRASIL into the radio regime, essentially following 
the prescriptions given by Condon &; Yin (1990). Before 
adopting it we have reviewed all the possible sources of 
radio emission related to the star formation activity, with 
particular emphasis on the integrated properties of stellar 
populations. 

In agreement with previous studies, we have found that 
the fraction of the NT emission due to radio supernova 
remnants is about 6% of the total and that other dis- 
crete sources provide a negligible contribution. Almost all 
the NT emission thus come from diffuse electrons possi- 
bly injected into the interstellar medium by CCSN events, 
and adopting a relation between the NT emission and the 
CCSN rate seems the safest way to proceed. However, con- 
trary to previous studies, we have also shown that the 



latter relation, which is at the base of the tightness of 
the FIR/Radio correlation, is a natural outcome whenever 
synchrotron losses happen on timescales that are shorter 
than the jading time of the CCSN rate. 

This is certainly the case in normal spirals, where the 
global SFR changes very little over the last billion years. 
But the situation may be very different in luminous star- 
bursts, where the SN rate may change significantly in a 
timescale typical of the lifetime of the most massive stars. 
Indeed, previous studies claimed that the existence of a 
FIR/Radio correlation under such conditions requires a 
fine tuning between the magnetic and radiative energy 
density, which is difficult to explain. We have shown that 
this fine tuning is not necessary. In fact the observed 
prevalence of synchrotron energy losses on inverse Comp- 
ton against the photons of the high stellar radiation field, 
indicates large magnetic fields and, consequently, guar- 
antees very short synchrotron electron lifetimes (Condon 
1992), certainly shorter than the typical lifetime of the 
most massive stars. The proportionality between the NT 
emission and the CCSN rate thus holds even in the ex- 
treme conditions found in the luminous obscured compact 
starbursts. This is why the FIR/Radio correlation appears 
so robust. 

We have calibrated the NT emission/CCSN rate rela- 
tion on the observed properties of our Galaxy. With these 
assumptions we reproduce well the FIR/Radio correlation 
of normal star forming galaxies, namely qi.4G# z =2.3. We 
are thus quite confident that our model is able to repro- 
duce both the FIR and Radio emission of star forming 
galaxies, with a minimum set of well calibrated parame- 
ters. 

As for the starburst galaxies, with a rapidly changing 
SFR, we have shown that the different fading times of the 
FIR and Radio emission may be used to analyse in great 
detail the recent star formation history in these galax- 
ies. Indeed, in the cases of M82 and ARP220, presented 
in Sect. |[ the combination of FIR and Radio observa- 
tions allows to reach a time resolution of a few tens of 
Myr, which is not possible based on the UV-FIR contin- 
uum properties alone, and even on the optical/NIR emis- 
sion lines diagnostics, for obscured galaxies (Poggianti et 
al. 2001, Bressan et al. 2001). We have thus analysed in 
greater detail the evolution of obscured starburst galaxies 
under different assumptions concerning the burst charac- 
teristics, challenged by the possibility of using the devia- 
tions from the FIR/Radio correlation as a diagnostic tool 
to infer the recent star formation history. We have com- 
pared our results with observations of a sample of compact 
ULIRGs, having in mind that, in these objects, it is also 
not clear what fraction of IR and radio emission is pos- 
sibly contributed by the AGN. Compact ULIRGs show a 
noticeable deviation from the average FIR/Radio relation, 
with their qi. 49GHz being generally lager than 2.3, a few of 
them with values as high as 3. Taken at face values, these 
deviations suggest that radio emission is depressed by a 
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large factor, relatively to normal spirals. Understanding 
the interplay between FIR and radio emission in these ob- 
jects is thus fundamental to make reliable predictions for 
high redshift dust enshrouded galaxies. 

Starburst models with peak SFR reaching several hun- 
dred M Q /yr and thereafter exponentially declining, may 
account for the IR and radio emission of the observed 
ULIRGs and are able to reproduce the observed variation 
of the value of qi. 49GHz • This view is consistent with the 
current idea that ULIRGs and, to a larger extent high red- 
shift dust enshrouded galaxies, are transient phenomena 
that nevertheless build up a significant fraction of stars 
and metals (Granato et al. 2001). 

The introduction of a new diagnostic diagram, the 
qi.49GHz vs radio slope diagram, allows us to single out 
the effects of starburst evolution and free- free absorption. 
Very young star bursts display an excess of FIR emis- 
sion relative to the radio emission because the latter is 
initially contributed mainly by the free-free emission pro- 
cess. As the starburst ages, the NT contribution increases 
and becomes the dominant source, while the radio slope 
reaches the typical values observed in synchrotron emis- 
sion. Free-free absorption affects the 1.49 GHz data, in- 
troducing a trend with the higher q being accompanied 
by the shallower slope. The estimated optical depths for 
free-free absorption at 1.49 GHz are between 0.5 and 1. At 
8.4 GHz, free-free absorption becomes negligible and the 
above trend disappears. The value of qs.4GHz is a measure 
of the age of the starburst. However, even in the latter 
diagram the slope is still affected by free-free absorption. 
Thus we suggest that a similar diagram between 8.4 GHz 
and a higher frequency range would be critical for the un- 
derstanding of the evolutionary status of compact ULIRGs 
because, in that case, the slope, unaffected by free-free 
absorption, would provide an independent estimate of the 
age. 

If ULIRGs are transient phenomena as suggested by 
other independent studies, then determining their SFR 
from conventional estimators may be a problem. They are 
far from being in a stationary status, the term "average 
star formation" is meaningless, and applying standard cal- 
ibrations may result in a significant error and/or appar- 
ent discrepancy between the observable themselves. One 
should be able to reconstruct the recent history of star 
formation and, for this purpose, we suggest the use of the 
above diagram to determine the characteristic parameters 
of the burst first, and then the age-averaged SFR from 
either the FIR or radio luminosity. 

Another relevant question addressed is how reliable is 
the use the FIR/Radio correlation to evaluate the con- 
tribution of non thermal radiation from the central ac- 
tive nucleus. Among the plotted data, the symbol "M" 
indicates the position of the Seyfert 1 galaxy Mrk 231 
(UGC 08058) (Thean et al. 2000). The fact that Mrk 231 
is clearly distinct from the other objects and occupies a 
position below any starburst model, becomes particularly 



evident in the qs.4GHz vs slope diagram, where the effects 
of free-free absorption on the q ratio are minimized. 

We have also shown that during the post starburst 
phase, the models reach values of q significantly lower than 
those of quiescent spirals, with still significant FIR lumi- 
nosities. This is consistent with the detection in nearby 
Abell clusters of a statistically significant excess of star 
forming galaxies with enhanced radio emission relative to 
the FIR (Miller & Owen 2001). We suggest that these low 
values of q are due to an evolutionary effect rather than a 
direct enhancement of radio emission by interaction with 
the intracluster medium. 

Finally we have investigated on the redshift depen- 
dence of the the FIR/Radio correlation and on its valid- 
ity (through the sub-mm radio index, s^ ) to provide a 
photometric redshift estimate of obscured distant galaxies. 
The large dispersion of sf 5 f observed among local galax- 
ies is compatible with the evolutionary effects discussed in 
Sect. ^|. The unknown evolutionary status of the starburst 
renders the sf 5 ^ index very unreliable at almost any red- 
shift. We thus suggest to complement the index sf 5 4 with 
a radio slope determination, because of its tight relation 
with the evolutionary phase of the starburst and the its 
very shallow dependence on the redshift. Other uncertain- 
ties like the sub-mm slope and the presence of a significant 
contribution at radio wavelengths from a central AGN, ob- 
viously worsen the above picture. 
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